The role of electronic correlation in the Si(100) reconstruction: 
a quantum Monte Carlo study 
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Recent low-temperature scanning tunneling experiments have challenged the generally accepted 
picture of buckled silicon dimers as the ground state reconstruction of the Si(100) surface. Together 
with the symmetric dimer model of the surface suggested by quantum chemistry calculations on 
small clusters, these findings question our general understanding of electronic correlations at sur- 
faces and its proper description within density functional theory. We present quantum Monte Carlo 
calculations on large cluster models of the symmetric and buckled surface, and conclude that buck- 
ling remains energetically more favorable even when the present-day best treatment of electronic 
correlation is employed. 
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Despite extensive experimental and theoretical inves- 
tigation, the nature of the reconstruction of the Si(100) 
surface is still subject to debate. While this surface is 
of technological relevance because of its use in the fabri- 
cation of electronic devices, determining its ground state 
structure is important as a test of our general under- 
standing of the role of electronic correlation at surfaces. 

Scanning tunneling microscopy (STM) experiments 
indicate that Si(100) reconstructs in rows of silicon 
dimers 0. At room temperature, most dimers appear 
symmetric due to their dynamical flipping motion, and 
only dimers close to defects are pinned in a buckled con- 
figuration H . The number of symmetric dimers decreases 
below 120 K and dimers buckle alternately within each 
row, with the formation of p(2x2) or c(4x2) domains cor- 
responding to adjacent rows in identical or opposite ori- 
entations [0|§- Experimentally, the c(4 x 2) reconstruc- 
tion (Fig. P was accepted as the lowest energy structure. 




FIG. 1. Model of the silicon (100) surface in the c(4 x 2) 
reconstruction. 

This picture is now being challenged by new experi- 
mental work. A series of low temperature STM studies 
have recently reported that, while the c(4 x 2) structure 
is observed below 120 K, further cooling below 20 K 
causes the dimers to appear again symmetric. Yokoyama 



and Takayanagi Q argue that this is a dynamical phe- 
nomenon caused by a lowering of the potential energy 
barrier between the two buckled configurations, which 
allows the dimers to resume the flip-flop motion charac- 
teristic of room temperature. In contrast, Kondo et al. || 
claim that the observed symmetric dimers are static since 
their images do not exhibit the noise associated with the 
flipping motion observed in the same sample at 110 K. 
However, evidence has also been produced that, at low 
temperatures, the rate of dimer flipping is significantly 
affected by the amount of tunneling current |6|], casting 
doubts on the conclusions drawn by Kondo. 

To date, theoretical work also remains divided on the 
issue. Density functional theory (DFT) calculations on 
slab geometries and large cluster models |5MITo|] fa- 




vor a buckled reconstruction. GW calculations 111] 
on buckled geometries show good agreement with the 
measured dispersion of surface band states JlH , and 
also a surface core-level shift analysis [Q supports the 
DFT finding. On the other hand, multiconfiguration 
self-consistent field (MCSCF) and configuration interac- 
tion (CI) Jl5| , p~6 17, 1^] calculations on small clusters find 
the symmetric reconstruction to have the lowest energy. 
These quantum chemistry techniques emphasize differ- 
ent components of electronic correlation than DFT: static 
correlation, arising from near-degeneracy of molecular or- 
bitals, is most effectively described by a linear combina- 
tion of low-lying determinants as in a MCSCF calcula- 
tion, while dynamical correlation, given by short-range 
electronic screening, is adequately treated in DFT. 

The buckling of the Si dimers on the Si(100) surface 
is one of the hard problems of many-body physics at 
surfaces, because it involves subtle aspects of electronic 
correlation JT9j ] . Since the symmetric surface dimer has 
two dangling bonds, each filled with one electron, it is in 
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a bi-radical state, and static correlation must be prop- 
erly taken into account. The buckled state, on the other 
hand, is thought to be stabilized by a rehybridization of 
the dangling orbitals, accompanied by charge transfer: 
the lower Si atom approaches sp 2 -like hybridization and 
its charge is depleted in favor of the higher Si atom in 
the dimer. Thus, the net stability of the buckled con- 
figuration depends on the degree to which the repulsive 
Coulomb interaction of the two electrons in the upper Si 
electronic state is dynamically screened. This is where 
dynamical correlation enters into the problem. 

In this Letter, we use quantum Monte Carlo (QMC) 
methods to find accurate energy differences between the 
symmetric and buckled reconstructions for large cluster 
models of the surface. Unlike the other theoretical meth- 
ods previously used for this problem, QMC has the ad- 
vantage that it can be applied to sufficiently large sys- 
tems and still provide an accurate description of both 
dynamical and static electronic correlation p(J . We find 
that dimer-dimer interactions are important and suffi- 
ciently large clusters must be used to adequately model 
the Si(100) surface. Our many-body calculations conclu- 
sively show that the ground state of the Si(100) surface is 
a buckled reconstruction, and that the trend with respect 
to cluster size found in DFT calculations is correct. 

Cluster models of Si(100). While methods based on 
DFT are known to usually give a good description of 
structural and elastic properties (e.g. surface lattice con- 
stant and elastic interaction between the Si dimers), it is 
unclear whether they can adequately represent the sub- 
tle aspects of electronic correlation at the Si(100) sur- 
face. It is therefore appropriate to use QMC on geome- 
tries obtained from DFT calculations to assess whether 
a more accurate treatment of electronic correlation can 
fundamentally change the picture. Here, we choose to 
address this issue by performing calculations for clusters 
that mimic the surface geometry. 

In identifying appropriate cluster shapes for Si(IOO), 
one must recognize that interactions are negligible be- 
tween neighboring dimer rows, while they are substan- 
tial between dimers in the same row. This is apparent 
from the small and the large dispersion of the surface 
band states along the respective directions, T — J and 
T—J' |ia] . Moreover, the p(2 x 2) and c(4 x 2) reconstruc- 
tions are found energetically quite close in experiments 
and calculations pl| (within 2 meV), demonstrating that 
dimer rows are weakly interacting. 

Therefore, the surface can be modeled with clusters 
containing only a single row of dimers. Such clusters 
with one, two and three dimers are SigHi2, SiisHig and 
Si2iH 2 o, previously also used in Refs. 
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They represent a four layer cut of the Si(100) surface with 
all but the surface atoms terminated with hydrogens to 
passivate dangling bonds. In Fig. 0, we show the Sii5Hi 6 
and Si2iH2o clusters both in the buckled and symmetric 
reconstructions. 



Most CI calculations have been performed on the 
smallest cluster SigHi2 which is found to be symmet- 
ric jHQ. Also for the Sii 5 Hi 6 cluster, Paulus Jl|] 
claims a symmetric ground state but the level of the CI 
calculation is not specified and the energetics are not ex- 
plicitly given. 

The SigHi2 cluster is however too small to draw con- 
clusions on the real surface, as became clear from recent 
DFT calculations ||. For SigHi2, approximate density 
functionals yield conclusions similar to CI, in that buck- 
ling is cither energetically unfavorable or marginally pre- 
ferred by less than 4 meV. However, as the number of 
dimers in the cluster increases, buckling becomes favor- 
able and the optimal buckling angle increases for all the 
functionals. For the three-dimer cluster Si2iH 2 o, the en- 
ergy gain per dimer is between 0.15-0.20 eV, depending 
on the functional used, and in agreement with the slab 
results to better than 0.05 eV. This suggests that it is 
possible to infer the behavior of the Si(100) surface from 
the three-dimer cluster. 




FIG. 2. S115H16 (upper figure) and Si2iH2o (lower figure) 
cluster models of the Si(100) surface. Both the symmetric 
(left) and buckled (right) reconstructions are shown. The 
bigger atom is silicon and the smaller one is hydrogen. 

Here, we perform QMC calculations on the clusters 
with two and three dimers, Sii5Hi6 and Si2iH 2 g. For 
both clusters, we use the geometries optimized within 
DFT using the PW91 functional |H|. For details on the 
construction and the geometry of the clusters, see Ref. ||. 

Silicon pseudopotential. For the silicon atom, we use a 
norm-conserving sp-non-local pseudopotential for the ten 
core electrons. The pseudo-potential was generated in 
an all-electron Hartree-Fock calculation for the Si atom. 
It was tested within QMC to calculate binding energy 
and bond length of Si2 which were found in excellent 
agreement with experiment. The transferability of the 
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pseudopotential with respect to all electron calculations 
was checked with Hartree-Fock (HF) and B3LYP || for 
larger silicon clusters. 

Quantum Monte Carlo methods. The form of the 
many-body wave function used in QMC calculations effi- 
ciently describes the static part of electronic correlation 
by the use of a linear combination of Slater determinants, 
as well as its dynamical component by introducing a pos- 
itive Jastrow correlation factor (modified from Ref. p4| ] 
to deal with pseudo-atoms): 

* = dn °n °i n j (»■« ' r » ' *» ■ ( x ) 

n aij 

Djj and are Slater determinants of single particle or- 
bitals for the up and down electrons, respectively, and 
the orbitals are represented using atomic Gaussian ba- 
sis |2f| ■ The Jastrow factor correlates pairs of electrons % 
and j with each other, and with every nucleus a, and dif- 
ferent Jastrow factors are used to describe the correlation 
with a hydrogen and a silicon atom. 

The determinantal part of the wave function is gener- 
ated within HF or MCSCF, using the quantum chemistry 
package GAMESS @. As active orbitals in the MCSCF, 
we choose the occupied bonding and the unoccupied an- 
tibonding ir orbitals, which are shown for the SiisHig 
symmetric cluster in Fig. 0. 




is tested using variational Monte Carlo (VMC). Reopti- 
mizing the determinantal part of the wave function pg| ] 
in the presence of the Jastrow factor did not lead to a 
significant improvement in the energy. The wave func- 
tion is then used in diffusion Monte Carlo (DMC), which 
produces the best energy within the fixed-node approxi- 
mation (i.e. the lowest-energy state with the same nodes 
as the trial wave function) p9| |. 

Results and conclusions. Since the interplay between 
static and dynamical correlations is central to this prob- 
lem, we first investigate the importance of accounting 
for near-degenerate molecular orbitals in the two-dimcr 
clusters. For the symmetric geometry, the inclusion 
of determinants corresponding to the active 7r-orbitals 
(Fig. ||) yields a better wave function, with a VMC en- 
ergy which is 0.34±0.06 eV per dimer lower than the 
single-determinant energy. However, within DMC, sin- 
gle and multi-determinant wave functions yield energies 
which only differ by 0.04±0.02 eV per dimer. Since the 
DMC energy of the symmetric reconstruction is rather 
insensitive to the use of more than one determinant, it 
is not surprising that, for the buckled geometry, no en- 
ergy gain is obtained neither in VMC nor DMC by using 
a multi-determinant wave function: the larger HOMO- 
LUMO gap of the buckled cluster makes a single deter- 
minant the best option also at the variational level. 

Since a multi-determinant treatment fails to yield a 
significant energy gain for the two-dimer cluster, the cal- 
culations for the Si2iH2o cluster are performed with one 
determinant only. In Table |, we list the QMC results 
obtained for the two- and three-dimer clusters using a 
single-determinant wave function. 

TABLE I. VMC and DMC energy differences between the 
symmetric (-B sy m) and buckled (£Wck) reconstructions of the 
S115H16 and Si2iH2o clusters. AE — E syln — -Ebuck and the 
numbers in parenthesis are the statistical errors on the last 
two figures. Energies are in eV. 







AE 


Ai?/dimer 


SilsHie 


VMC 


0.53(13) 


0.27(06) 




DMC 


0.06(04) 


0.03(02) 


Si2iH2o 


VMC 


0.70(14) 


0.23(05) 




DMC 


0.34(06) 


0.11(02) 



FIG. 3. Bonding (upper figure) and antibonding (lower fig- 
ure) 7r orbitals for the symmetric SiisHi6. The orbitals can 
be bonding (right) or antibonding (left) between adjacent 
dimers. 

Once the determinantal part of the wave function has 
been determined, the parameters in the Jastrow factor 
are optimized within QMC using the variance minimiza- 
tion method p7j] and the accuracy of the wave function 



For Sii5Hi6, the symmetric and buckled reconstruc- 
tions are energetically very close as in the case of the 
one-dimer cluster. With a DMC energy difference per 
dimer of 0.03±0.02 eV, we cannot establish which recon- 
struction is more favorable for this cluster size. However, 
for Si2iH 2 o, the DMC energy gain per dimer in favor of 
buckling increases substantially to 0.11±0.02 eV. This 
is in good agreement with the trends established within 
DFT that a one- or a two-dimer cluster fails to model the 
surface accurately M. 
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In Fig. |, we summarize the energy differences between 
the symmetric and the buckled reconstruction obtained 
within the local-density approximation (LDA), several 
generalized gradient approximations (GGA), B3LYP and 
DMC. We also include results for SigHi2 and the slab 
geometry ||. DFT functionals predict that the buck- 
led cluster is favored by 0.07-0.13 and 0.15-0.20 eV per 
dimer for Sii5Hi6 and Si2iH2o, respectively. Even though 
the DMC energy differences are noticeably smaller than 
in DFT, they clearly indicate that buckling is energeti- 
cally more favorable. The DMC results also highlight the 
importance of dimer-dimer interactions: SiisHig behaves 
still very much like the single-dimer cluster and a Si2iH2o 
cluster is needed to be able to resolve the energy gain in 
favor of the buckled geometry. 
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FIG. 4. Energy differences per dimer for the symmetric 



and buckled reconstructions (AE = E s , 



£Wck) obtained 



with LDA, various GGA functionals, and DMC. 

In this Letter, we presented accurate QMC calculations 
of the energies of the buckled and symmetric reconstruc- 
tion of the Si(100) surface. Using large clusters, we 
determined that the buckled geometry is lower in energy. 
Thus, recent speculations fueled by low-temperature 
STM experiments that electron correlations at low tem- 
peratures would favor symmetric dimers can be ruled 
out. Furthermore, our calculations show that electronic 
interactions between adjacent Si dimers in a row are im- 
portant to obtain the buckled ground state, and large 
clusters must be used to adequately model the Si(100) 
surface. While DFT-LDA/GGA calculations tend to 
overestimate the energy gain due to buckling, the cor- 
rect trend with cluster size is reproduced already at the 
LDA/GGA level of theory. 
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